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PARAMETER  ESTIMATION  TECHNIQUES 
for 

NONLINEAR  DISTRIBUTED  PARAMETER  SYSTEMS 
H.T.  Banks  and  K.  Kunisch 


ABSTRACT 

Methods  for  estimating  system  parameters  are  discussed  for 
a  class  of  partial  differential  equations.  We  develop  schemes 
based  on  modal  subspace  approximations  in  some  detail  and  include 
numerical  examples. 
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1 .  Introduction 

In  this  paper  we  study  approximation  methods  for  linear  and 
nonlinear  partial  differential  equations  and  associated  parameter 
identification  problems.  As  will  be  seen,  many  of  the  ideas  are 
classical  in  nature,  but  the  proposed  .operator-theoretic  approach 
is  appealing  for  its  conciseness  and  generality.  The  parameter 
identification  problem  is  explained  in  Section  2.  In  Section  3 
we  develop  general  nonlinear  approximation  results  which  are  sub¬ 
sequently  used  for  modal  approximations  of  certain  classes  of 
hyperbolic  and  parabolic  partial  differential  equations.  In  the 
fourth  section  the  theory  developed  so  far  is  employed  to  approxi¬ 
mate  the  (infinite-dimensional  state)  parameter  identification 
problem  by  a  sequence  of  problems  for  ordinary  differential 
equations.  We  have  tested  the  modal  approximation  scheme  in 
numerical  experiments  for  parameter  identification  in  hyperbolic 
equations;  some  of  our  findings  are  reported  in  Section  5.  De¬ 
tailed  proofs,  further  numerical  results,  and  a  discussion  of 
relevant  literature  on  these  problems  will  appear  in  a  forthcoming 


paper. 


2. 


2 .  The  Identification  Problem 

We  consider  the  abstract  Cauchy  problem 

u(t)  -  A(q)u(t)  +  F  (q  ,t ,u(t) )  ,  t  >  0 

(2.1) 

u(0)  *  uQ(q) 

where,  for  each  q  €  Q  c  R  ,  A(q)  is  the  infinitesimal  generator 
of  a  linear  Cg-semigroup  (T(t;q)}t>0  on  a  Hilbert  space  X(q). 
The  inner  product  and  norm  in  X(q)  will  be  denoted  by  (•,•) 
and  |*|  respectively,  although  on  occasion  we  shall  omit  the 
subscript  q.  Recall  that  t  -*•  u(t)  is  called  a  mild  solution 
of  (2.1)  if 


u(t) 


T(t;q)u0 


T(t-o;q)F(q,o,u(a))do. 


Conditions  on  F  that  guarantee  existence  of  mild  solutions  to 
(2.1)  will  be  given  below.  For  the  relationship  between  mild  and 
strong  solutions  of  (2.1)  one  may  consult  [4  ];  here  we  only  note 
that  in  many  specific  instances  (i.e.,  hypotheses  on  A(q)  and/or 
F  and/or  Uq)  mild  solutions  are  in  fact  strong  solutions. 

We  assume  in  our  discussions  that  X(q)  is  a  function  space 
of  Rn-valued  "functions"  (or  the  usual  Lebesgue  equivalence 
classes)  defined  on  a  fixed  interval  [0,1]  and  thus  we  shall 
also  use  the  notation  u(t,x;q)  or  u(t,*;q)  to  denote  solutions 
of  (2.1).  We  shall  present  approximation  techniques  for  para¬ 
meter  identification  problems;  these  will  be  discussed  in  the 
context  of  a  typical  least  squares  problem.  For  example,  at 
points  (Xj),  0  <_  x^  <  Xj  <  *•*  <  x^  <_  1,  and  times  {t^}, 

0  -  fci  <  t2-“  <  tr  <  T,  we  might  be  given  s-vector 
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measurements,  s  £  n,  of  the  n-dimensional  "state*',  which  are 
denoted  by  y\  €  R^,  i  =  l,...,r.  These  are  assumed  to  represent 
measurements  of  C(q)£(t^;q)  where  5(t^;q)  =  col (u(.t^ ,x^;q)  , . . . , 
u(ti,x^;q))  and  C(q)  is  an  (sfc)  x  (ni,)  matrix  depending  con¬ 
tinuously  on  q.  A  typical  identification  problem  (ID)  is  then  given  by: 

(ID) :  Minimize 

J(q)  =  .1  |C(q)C(ti;q)-yi|2 

over  q  €  Q. 

We  note  that  often  one  cannot  observe  all  components  of  the  state 
u  and  hence  it  is  necessary  (in  order  to  conform  with  reality) 
to  use  the  matrix  C(q)  in  the  problem  formulation  and  analysis. 

We  also  remark  that  the  point  evaluations  (at  Xj)  used  to  define 
Cftpq)  above  may  be  meaningful  only  in  specific  instances,  depend¬ 
ing  on  the  equation  (2.1),  the  space  X  and  the  initial  data  Ug. 

For  the  special  cases  presented  below,  one  can  verify  that  the  mild 
solutions  we  consider  do  yield  functions  for  which  point  evaluation 
is  a  valid  operation.  The  least  squares  criterion  that  is  used  in 
defining  (ID)  is  just  one  of  several  possible  choices  of  criteria 
that  might  be  used  in  identification  or  parameter  estimation  problems 
and  it  will  be  obvious  that  our  discussions  and  analysis  here  can 

be  extended  to  cover  other  types  of  functionals  J  which  satisfy  hypothesis 
(H10)  of  vSection  4. 

As  indicated  above  we  approximate  (ID)  by  a  sequence  of  problems 
(ID^),  each  of  which  can  be  solved  with  standard  numerical  procedures. 


4. 


To  formulate  the  approximating  problems,  we  take  for  each  q  a 

M 

sequence  of  closed  linear  subspaces  X  (q)  of  X(q)  endowed  with 

the  topology  induced  by  X(q).  The  orthogonal  projections  of 

X(q)  onto  X^(q)  are  denoted  by  ir^(q).  We  then  define  the 

operators  A^(q):X(q)  -*  XN(q)  approximating  A(q)  by 

AN(q)  =  TiN(q)A(q)TrN(q)  where  Dom(A^(q))  »  X(q) .  This  form  of 

approximating  operators,  which  is  a  classical  one  (e.g.,  see 

[  5  ,  p.369]),  has  recently  proved  to  be  useful  in  deriving  approximation 

methods  for  delay-differential  equations  [1  ].  We  stress  that  this 

formulation  entails  the  implicit  assumption  XN(q)  <=  Dom(A(q)). 

The  projections  of  F  onto  the  subspaces  are  given  by 
FN(q,t,v)  =  ir^(q)F(q,t  ,v)  for  each  (q,t,v)  €  Q*[0,T]xX  (we  shall 
assume  throughout  that  the  spaces  X(q)  are  set-wise  all  the  same 
set  X).  The  family  of  approximating  equations  is  therefore 
given  by 

v(t)  =  AN(q)v(t)  +  FN(q,t,v(t)) ,  t  >  0 

(2.2)  N 

v(0)  “  TiN(q)u0(q). 

The  hypotheses  on  F  given  below  will  insure  existence  of 
mild  solutions  u^(the  notation  u^(t) ,u^(t ;q)  and  uN(t,x;q) 
will  all  be  used  in  the  sequel)  of  (2.2).  Since  X^(q)  is  in¬ 
variant  under  AN(q),  (2.2)  is  easily  seen  to  be  an  initial  value 
problem  in  the  subspace  XN(q).  In  the  event  that  X^(q)  is 
finite-dimensional,  (2.2)  is  equivalant  to  a  system  of  ordinary 
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differential  equations  for  the  generalized  Fourier  coefficients 
of  the  representation  for  v  relative  to  a  chosen  basis  for 
XN(q).  In  our  discussions  we  make  the  assumption: 

(HI)  All  elements  of  XN(q)  are  piecewise  continuous  functions 
on  [0,1]. 

Defining  CN(ti;q)  =  col (uN(ti ,x1 ;q) . uN(ti ,xA;q) ) ,  we 

formulate  the  approximate  identification  problems  (ID^)  corresponding 
to  (ID)  by: 

(ID1^)  :  Minimize 

jN(q)  =  Ji|C(qHN(ti;q)-yi|2 

over  q  €  Q. 

Before  discussing  the  existence  of  solutions  to  problems  such  as  (ID)  and  (ID^) 
and  their  relationship,  we  present  convergence  results  for  the 
approximation  of  (2.1)  by  (2.2)  in  a  form  readily  applicable  to 
the  identification  problems. 

We  shall  call  (q  ,XN(q)  ,irN(q)  ,  AN(q)  ,FN(q)  }  an  approximation 

scheme  for  (2.1)  if  q  €  Q  c  Rk,  XN(q)  is  a  sequence  of  subspaces 

of  X(q)  and  irN(q)  ,AN(q)  ,FN(q)  are  maps  *N(q):X(q)  -*-XN(q), 

AN(q):X(q)  -»■  XN(q)^N(q)  :  [0  ,T]  *  X(q)  -*■  XN(q).  Such  a  scheme  will 

be  said  to  be  convergent  if  for  any  qN,q  6  Q,  lim  qN  *  q  implies 

N-*-» 

that  corresponding  mild  solutions  u^(t;q^)  of  (2.2)  converge  to 
a  mild  solution  u(t;q)  of  (2.1). 


3.  Modal  Approximations:  Convergence  Results 

In  this  section  we  present  a  convergence  theorem  for  non¬ 
linear  systems  (2.1),  (2.2)  and  then  discuss  modal  approximation 
schemes  for  classes  of  hyperbolic  and  parabolic  equations.  We 
shall  refer  to  a  number  of  hypotheses  which  we  present  now. 

(H2)  For  each  q  €  Q,  A(q)  generates  a  linear  (^-semigroup 
T(t;q) . 

(H3)  The  spaces  X(q),  q  €  Q,  are  set-theoretically  equal  and 
AN(q)  generates  a  C  -  semigroup  TN(t;q)  on  X(q). 
Moreover,  q^  -*■  q  implies  |T^(t  ;q^)z  -  T(t;q)z|  ^  0 

q 

as  N  ®  for  all  z  €  X(q),  uniformly  in  t  on  compact 
subsets  of  [0,“).  There  exist  constants  M  and  to  independent 
of  N  and  q  such  that  |T(t;q)|  £  Meut  and  |TN(t;q)|  <_  Mewt. 

(H4)  The  set  Q  <=  R  is  compact. 

(H5)  The  spaces  X(q)  are  uniformly  topologically  isomorphic. 
That  is,  there  exists  a  real  constant  K  such  that 
|vU  <  K | v | q  for  all  q,q  in  Q. 


(H6)  The  projections  irN(q):X(q)  XN(q)  satisfy:  For  any 
sequence  {q^>  in  Q  with  q^  q,  one  has 
|ir^(qN)z  -  z  |  m  0  as  N-*®  for  each  z  6  X(q). 

q 
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(H7)  The  nonlinear  function  F : Q  x  [0,T]  x  X(q)  -*•  X(q)  satisfies: 

(i)  For  each  continuous  function  u:[0,T]  -*■  X(q), 
the  map  t  -*•  F(q,t,u(t))  is  measurable  for 
each  q  e  Q. 

(ii)  There  exists  a  function  kj  in  L2(0,T;R)  such 

that 

|F(q,t,u^)  -  F  (q,t  ,U2)  |  q  5  ^i  (t)  |  u-^-u2  I  q 
for  all  q  €  Q,  u1,u2  €  X(q). 

(iii)  There  exists  a  function  k2  in  L2(0,T;R)  such 

that 

|F(q,t,0)|q  <  k2(t) 

for  all  q  €  Q. 

(iv)  For  each  (t,u)  €  [0,T]  x  X(q),  the  map 
q  F(q,t,u)  is  continuous. 

We  remark  that  to  be  more  precise,  we  should  have  written 
|ir^(q^)z  -^Nz|  jyj  0  *  where  is  the  canonical  isomorphism 

—  Q  XT 

of  X(q)  onto  X(q  ),  in  (H6)  (a  similar  adjustment  should  be 
made  in  (H3) ).  However ,  we  suppress  this  notation  throughout  in 
as  much  as  there  will  be  no  confusion  in  light  of  the  assumed  set¬ 
wise  equality  of  the  X(q),  q  €  Q.  We  further  remark  that  (H5) 
implies  uniform  (in  q  and  N)  boundedness  of  |ir^(q)|  . 

Theorem  3.1:  Suppose  that  (H2)  -  (H7)  obtain  and  that  {qN}  is_ 

N  — 

any  sequence  in  Q  with  q  ■+■  q*  Further  assume  that 
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|uQ(qN)  -  Ug(q)|  N  -*■  0  a£  N  °°.  Then  for  each  q  €  Q,  unique 

M 

mild  solutions  u(t;q)  and  u  (t;q)  of  (2.1)  and  (2.2)  exist 
and  we  have  |u^(t;q^)  -  u(t;q)|  N  0  a£  N  -*•  °°  for  each 

t  e  [o ,t]  .  q 

Corollary  5.1:  If  (H7)  is  strengthened  so  that  one  assumes 
(o,v)  ■*  F(q,o,v)  is  continuous  on  [0,T]  x  X(q),  then  the  con¬ 
vergence  u^(t;q^)  -*■  u(t;q)  of  Theorem  3.1  is  uniform  in  t  on 
[0,T]. 

Example  3.1.  (Hyperbolic  equations). 

We  consider  here  the  equation 

(3.1)  utt  =  qxuxx  +  q2ut  +  q3u  +  f  (qg . qJJ.t ,  •  ,u) 

with  initial  and  boundary  conditions 

m  . 

u(0,x)  *  I  qj^i(x), 

(IC) 

m  . 

«t(0,x)  =  l  qr<l'i(x),  0  <  x  <  1, 
t  i-i  b  1 


(BC) 

u(t,0)  =  u(t,l)  =  0, 

t  >  0, 

where 

u(t,x)  and  q^,qj  are  scalars. 

9. 


Remark  5.1:  Although  in  (3.1),  (IC),  (BC)  we  do  not  explicitly 
allow  nontrivial  boundary  conditions  (possibly  containing  parameters), 
it  can  easily  be  seen  that  such  situations  are  included  in  our  for¬ 
mulation.  For  consider 

(3-2)  utt  "  Vxx 


with  the 

(IC) 

(BC) 


initial  conditions  and  boundary  conditions 

u(0,x)  =  q4<t>(x) 

«t(0,x)  =  q5iKx) 

u(t  ,0)  =  qyb^t),  u  (t ,  1 )  =  qgb2(t), 


where  bj,b2  aTe  twice  continuously  differentiable  functions. 
Then  the  usual  transformation  given  by  w(t,x)  *  u(t,x)  - 
(l-x)q^bj(t)  -  xqgb2(t)  transforms  (3.2),  (IC),  (BC)  into  the 
problem 


Wtt  *  «lwxx  -  (1-It)l7bltt  -  xc*8b2tt 
w(0,x)  =  q4<Kx)  -  (l-x)q7b1(0)  -  x  qgb2(0) 

wt(0,x)  =  q5iHx)  -  (l-x)qyblt(0)  -  xq^b2t(0) 
w(t  ,0)  =  w(t,l)  =  0, 


which  is  a  problem  that  is  a  special  case  of  the  formulation  (3.1) 
(IC),  (BC)  above. 

To  treat  (3.1),  (IC),  (BC)  we  first  rewrite  (3.1)  as  an  ab¬ 


stract  equation  in  the  usual  manner  employing  the  operator  A  = 


_ 


3x 


in  =  I^COjljR).  The  Sobolev  spaces  H1  considered  here  and 

below  will  consist  of  R  -valued  functions  on  [0,1]  and  we  suppress 

1  2 

this  notation  hereafter.  With  A  defined  onv-Dom(A)  =  Hq  n  H  ,  we 
include  the  boundary  conditions  (BC)  and  are  also  able  to  establish 
that  A  is  self-adjoint  and  (-Au,u)  |u|  for  all  u  €  Dom(A). 

Standard  results  guarantee  existence  of  A*^2  on  Dom(A^2)  =  Hq. 

We  make  the  following  additional  assumption: 


a  b  k 

(HQ)  There  exist  positive  numbers  q^,q^  such  that  q  €  Q  c  R 
implies  qa  ±  q1  < 


Then  the  set  Dom  (A1^2)endowed  with  the  inner  product  (u,w)  = 

ql 

(qlux’wx^L  a  Hilbert  space  which  we  donote  by  V(q^).  Finally, 

H(q^)  =  V(q^)  *  H^  with  the  product  topology  is  a  Hilbert  space 
X(q)  in  which  we  can  now  rewrite  (3.1),  (IC),  (BC)  as 


d  /uCt)\  /«(t)v 

jF  (  )  -j/t q)(  )  ♦  F(q,t,u(t)),  t  >  0 

dt  Vv(t);  Vv(ty 


where  (4k, ^i)  €  Hfqj) ,  Dom(.o^(q))  *  (hJ  n  H2)  *  hJ, 


( 


,u(t ,  •)) 


/  0  1  \ 

Mq)  *  (  )  ,  and  F(q,t  ,u(t)) 

V»1A  +  q3  q2/ 


11. 


Here  q^  =  (qj,...,qj)  for  j  =  4,5,6  and  q  =  (qx ,q2 .q3  .q4 »q5 »q6) 

^TTI  4*  7 

is  restricted  to  the  fixed  set  Q  c  RJ 

We  see  immediately  that  in  view  of  assumption  (HQ),  the  spaces 
H(q^)  are  all  topologically  isomorphic;  indeed  for  any  pair  q^.q^ 
in  [q^qj]  we  find  |z|H£~^  <  K|  ^  I  H(q1)  with  K  =  Cq^/q^) 1/2 
so  that  (H5)  is  satisfied.  We  discuss  (H3) ,  (H6)  for  a  specific 
choice  (so-called  modal  approximations)  of  finite  dimensional  sub¬ 
spaces  and  operator  approximations.  We  refer  the  reader  to  [  5, 
pp.  247-250]  for  the  relevant  background  material  needed  in  our 


development . 

The  operator  A  is  selfadjoint  with  compact  resolvent  and 

^  oo  00  J’JF 

the  eigenvectors  {4.}.  ,and  {$>. }  where  4>.  (x)  =  sin  jirx 

J  T1  J  J  J" 

and  4^  (x)  *  /I  sin  jirx  constitute  complete  orthonormal 

sets  (CONS)  for  V(l)  =  and  H®  respectively.  We  define 

the  modal  subspaces  X^(q)  «  H^(q^)  of  H(q^)  by 


HN(q1)  -  span  | (  ),...,(  ),(  ),...,(  )|  . 

1  10  V0  n*n'J 


We  note  that 


u{(\(°)j 

j  =  l(  0  9.  ) 


forms  a  CONS  for 


H(l)  and  a  complete  orthogonal  (but  not  normal)  set  for  H(q^), 

qx  t  1. 


The  modal  approximations  .Q^N(q)  for  ja^(q)  are  defined  by 
j*fN(q)  =  7rN(q1  Xoitq) Tr^Cq^ ) »  where  trN(q1)  is  the  canonical  projection 
of  H(q^)  onto  H^(q^).  The  requirement  that  H^(q^)  <=  Dom(jj<(q)) 
is  trivially  seen  to  be  true  here. 
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Theorem  3.2.  Let  (HQ)  obtain  and  let  q  »q  €  Q  e  R^m+^  be 

N  —  _  N  N 

such  that  q  -*■  q  a_s  N  <».  Then  _Qf(q)  and  (q  )  generate 

Cq- semigroups  T(t;q)  and  TN(t;qN)  on  H(q^)  and  H(q^)  res¬ 
pectively.  Furthermore ,  there  is  a  constant  w  €  independent 
of  N  such  that  |T(t;q)  |  <_  eu)t ,  |TN(t;q^)  |  <  eut  for  t  ^  0, 
and  for  each  z  €  H(q^)  ,  |TN(t;qN)z  -  T(t  ;q)z  |  0  a£  N 

q 

uniformly  in  t  on  compact  subsets  of  [0  ,<*>) . 

The  proof  of  this  theorem  can  be  given  employing  a  generalized 
version  of  the  Trotter-Kato  theorem  (see  [3  ])  and  the  spectral 
theorem.  One  can  readily  establish  dissipativeness  of  _a^(q)  -  ojI 
in  H(q)  for  anappropriately  chosen  w.  A  consequence  of  (H5)  is 
that  one  actually  obtains  T^(t;q^)z  •+■  T(t;q)z  in  H(l). 

Turning  next  to  the  nonlinear  equation  (3.3),  we  let 
Q  »  (q^  €  Rm | q  €  Q}  and  make  the  following  hypotheses  on  f. 

(H7*)  The  nonlinear  function  f:Q  x  [0,T]  x  [0,1]  x  R*  R*  satisfies 

(i)  For  each  (q^,u)  €  Q  x  r*,  the  map 
(t,x)  -*■  f(q^,t,x,u)  is  measurable. 

(ii)  There  exists  k^  in  L2(0,T;R)  such  that 

|f(q6,t,x,u1)-f(q6,t,x,u2) |  <  k1(t)|u1-u2| 
for  all  q6  €  Q,  t  €  [0  ,T] ,  x  €  [0 ,1] ,  and  Uj,u2  €  R1. 

(iii)  There  exists  k2  in  L2([0,T]  x  [0,1] ;R)  such  that 

lfCq6»t,x,0) |  <  k2(t,x) 
for  all  q^  €  Q. 

(iv)  For  each  (t,x,u)  in  [0,T]  x  [0,1]  x  R* ,  the  map 


q^  f(q^,t,x,u)  is  continuous. 
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It  is  not  difficult  to  see  that  if  one  defines  F  by  F(q,t,z)  = 
col(0,f  (q,t,  •  ,u(0))  for  z  =  col(u,v)  in  Hfq^,  the  conditions 
(H7*)  for  f  imply  (H7)  for  F.  Furthermore,  it  is  rather  easy 
to  establish  that  the  projections  ir^(qj)  are  in  fact  independent 
of  q^  so  that  completeness  of  the  ^j»^j  al°ng  with  (H5)  yield 
(H6).  Since  (H2)  and  (H3)  follow  from  Theorem  3.2,  we  may  apply 
Theorem  3.1  to  obtain  convergence  of  the  modal  approximation  scheme 
{q,HN(q1)  ,  ttN  (qx)  ,_^N(q)  ,FN(q)  }  for  hyperbolic  systems  (3.1),  (IC)  , 
(BC). 

Example  3.2.  (Parabolic  equations) 

For  our  second  class  of  examples  we  consider  parabolic  equations 

ut  =  ^(pux)x  +  q2u  +  f^q4 . 

(3-4^  m  . 

u(0,x)  -  l  qi<J>* (x) ,  0  <  x  <  1, 

i=l  *  1 

subject  to  the  boundary  conditions 

(3.5)  RjU(t,  ■)  =  0  for  j  =  1,2 . 

Here  we  assume  <(>.  €  ,  u(t,x)  6  R1 ,  and  q  ■  (q^  ^  >^3  >^4)  with 

q.  *  (q* , . . . , q1?)  , for  j  =3,4.  The  operators  Rj  defining  the  boundary 

2 

conditions  have  domain  H  and  are  given  by 

RjV  ■  ajiv(0)  +  ®j2v' (0)  +  °j3v(1)  +  aj4v' (1) 

for  v  €  .  We  make  the  following  assumptions  on  k,p  and  : 
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(H8)  The  functions  p,px  and  k  are  in  C(0,1;R)  with  k(x)  >  0 
and  p(x)  >  0  for  0  <  x  <  1. 


(H9)  The  matrix 


“ll  a12  a13  a14 

a21  a22  a23  a24 


has  rank  2  and  we  have 


P(0)  {aiia22”a12a21}  =  P  (!)  {al.^0t24 ~al  da23 ^  • 


13a24"a14a23 


We  again  rewrite  our  problem  as  an  abstract  Cauchy  problem 
and  to  this  end  we  define  the  generalized  Sturm-Liouville  operator 
_Q^(q)  in  H°  by  Dom(Ja?(q))  =  (<*>€  H2|R^<J>  =  0,  j  =  1,2}  and 
ja^(q)<))  =  k  *(q^p<J>x)x  +  <12^*  Then  (3.4),  (3.5)  can  be  written  as 


u(t)  *  srf (q)u(t)  +  F(q,t,u(t)),  t  >  0, 


(3.6) 


m 

u(0)  =  l  qi*. , 
i=l  5  1 


where  £  H  and  F(q,t,u(t))  =  £(q^,t, • ,u(t , •)) •  We  consider 

this  equation  and  the  operator  ja^(q)  in  with  inner  product 

r  1 

<u,v)  =  u(x)v(x)k (x)dx  and  note  that  (unlike  the  formulation  for 
J  0  n 
Example  3.1)  our  Hilbert  space  X(q)  *  H  is  independent  of  q 

in  this  case. 

The  equality  in  (H9)  implies  that  J3^(q)  is  selfadjoint  and 

spectral  results  for  _c/(q)  (e.g.  see  [  2])  yield  existence  of  a 

o  00  ~  „ 

CONS  (in  H  )  of  eigenfunctions  {¥.}  of  j#(q)  where 

jMl 

q  =  (1,0,..., 0)  6  R2m+2.  As  in  Example  3.1,  we  define  the  approximat¬ 
ing  modal  subspaces  of  H°  by  HN  ■  span  . and  let 
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ttN  :  H°  -*  HN  denote  the  canonical  projections.  This  determines 
in  the  same  manner  as  before  the  operators  (q)  *  tt  _Q^(q)it  and 
FN  *  ir^F.  Using  the  theory  of  general  Sturm- Liouville  operators 
and  the  Trotter- Kato  theorem  one  can  establish  convergence  of  the 
corresponding  semigroups. 


Theorem  3.3.: 

N  _ 

q  ,q  €  Q  c  R 


Suppose  (HQ) ,  (H4) ,  (H8)  and  (H9)  hold  and  let 
l2M+2  be  such  that  qN  -*•  q  as  N  -*■  «°.  Then  stf (q) 

_  •  /  M  _  1  fpN  f  4*  •  fl  ^  \ 


and  _c^N(qN)  generate  CQ- semigroups  T(t;q)  and  T  (t;q  )  on 


H°  that  satisfy  | T(t  ;q)  |  <  ewt ,  |TN(t  ;qN)  |  <  e10  for  some  w 
which  is  independent  of  N.  Furthermore  TN(t;qN)z  -*  T(t;q)z  for 
each  z  €  H°  with  the  convergence  uniform  in  t  on  compact  sub¬ 
sets  of  [0,~). 


If,  in  addition  to  the  hypotheses  of  Theorem  3.3,  one  assumes 
that  f  satisfies  (H7*)  (with  q6  replaced  by  q4) ,  then 
Theorem  3.1  holds  to  yield  convergence  of  the  approximating  solutions 
of  (3.6)  and  thus  the  modal  approximation  scheme  (q,HN,TTN,ja/N(q)  ,FN(q)  } 
is  a  convergent  scheme  for  (3.4),  (3.5). 


4 .  Approximation  of  the  Identification  Problem 

Returning  to  identification  problems  such  as  those  discussed  in  -Section 
2,  we  are  now  in  a  position  to  establish  convergence  of  solutions  of  the 
approximate  problems  ( )  of  minimizing  (q)  =  ^/’(uN  (q) , q)  to  those  of 
the  problem  (£?)  of  minimizing  J(q)  =  ^/(u(q),q)  where  we  make  the  follow¬ 
ing  assumption  on  the  fit-to-data  functional 

(HLO)  The  mapping  fx  X(q)  *Q  +  R  is  a  continuous  functional;  here  X(q) 
is  endeared  with  any  of  the  equivalent  topologies  hypothesized  in  the 
standing  assumption  (H5) . 


Our  results  are  stated  precisely  in  the  form  of  a  theorem. 


16. 


Theorem  4.1.  We  assume  the  hypotheses  (HI)  -  (H7)  and  (H10)  hold  and  that 
N  N 

q  -+■  uQ(q),q  -*•  tt  (q)z  and  q  T  (t;q)z  are  continuous  for  each 
z  e  X(q)  and  t  e  [0,T].  Then 

(a)  for  each  N  there  exists  a  solution  qN  of  ( &„)  , 

N,  N 

(b)  there  exists  a  subsequence  {q  J  }  of  {q^}  con¬ 
verging  to  some  q  €  Q  which  is  a  solution  of  (^). 

N.  N . 


Move over  t 
.N 


u  J(t;q  J)-u(t;q) 


Ni 
q  J 


as  Nj  0 


where  u  ,u  are  solutions  of  (2.2),  (2.1)  respectively. 


Proof :  It  is  not  difficult  to  argue  (using  (H3) - (H7) and  (H10))  that 

q  J^(q)  is  continuous  on  the  compact  set  Q.  Existence  of  a 
Ni  Ni  - 

subsequence  {q  J }  with  q  J  q  follows  from  (H4).  We  observe 

Ni  Ni  Ni 

that  for  any  q  €  Q,  one  has  J  ^ (q  J)  <  J  J (q) .  From  Theorem 

Ni 

3.1,  (H5)and  (H10)  ,we  see  that  J  J  (q)  -*■  J (q)  for  each  q  e  Q  and 
Ni  Ni  - 

J  J  (q  J)  J(q)  so  that  J(q)  <  J(q)  for  any  q  €  Q.  That  is, 


q  is  a  solution  of 


If  we  further  assume  (HQ),  it  is  quite  simple  to  see  that 
Theorem  4.1  is  applicable  to  identification  problems  for  Examples 
3.1  and  3.2.  We  have  conducted  numerical  investigations  for  the 
modal  approximation  scheme  for  identification  problems  with  the 
hyperbolic  systems  of  Example  3.1  and  report  briefly  on  some  of 
them  in  the  next  section.  For  these  calculations  we  chose 
C(q)  =  (1,0)  in  the  functionals  J  and  JN  of  (ID)  and  (IDN)  of 
.Section  2,  thereby  enabling  one  to  verify  (H10) . 
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5 .  Numerical  Results 

In  this  section  we  present  two  examples,  deferring  a 
thorough  discussion  of  our  numerical  experience  with  modal 
schemes  to  the  forthcoming  paper.  A  standard  IMSL-package 
employing  the  Levenberg-Marquardt  algorithm  was  used  to  solve 
the  parameter  identification  problems  for  the  approximating 
ordinary  differential  equations.  "Exact"  solutions  for  the 
distributed  systems  below  were  generated  independently  by  a 
simple  Crank-Nicolson  algorithm  and  some  of  these  values  were 
used  for  the  "data"  y^  in  (IDN) .  We  would  like  to  express 
our  deep  appreciation  to  James  Crowley  for  his  efforts  in  develop¬ 
ing  the  software  packages  employed  in  our  computational  experiments. 

Example  1 .  We  consider  equation  (3.1)  with  f  =  0: 

utt  *  ql“xx  +  q2ut  +  q3u. 

where  u(0,x)  =  q^  x(l-x),  ufc(0,x)  =  2qgx  for  0  <_  x  £  5-  and 
ufc(0,x)  =  q5(2-2x) ,  for  j  <  x  <  1,  and  u(t,0)  =  u(t,l)  =  0 
for  all  t.  The  true  model  parameters  were  taken  to  be 
q  =  (1.414,0,0,4,5) . 

N  0 

We  performed  a  five  dimensional  search  starting  at  q  '  = 

(1,0, 0,1,0)  for  each  value  of  N;  the  results  obtained  by  apply¬ 
ing  the  Levenberg-Marquardt  procedure  for  several  values  of  N 
are  given  in  Table  1. 


Table  1 


N 

N 

q2 

q3 

0.0026 

-0.2762 

-0.0035 

-0.0308 

-0.0001 

-0.0025 

0.0000 

-0.0006 

O 

• 

O 

0.0 

Example  5.2.  In  this  example  we  present  computations  for  the 
nonlinear  equation 

utt  -  1luxx  +  «2ut  +  «3U  +  ITu  ' 


with  u(0,x)  =  q^x(l-x) ,  u^(0,x)  =  q^  and  homogeneous  boundary 
conditions.  Table  2  depicts  the  results  corresponding  to  the 
"true"  model  parameters  q  =  (1.414,0,1,2,0).  Holding  q2 

and  q,.  fixed,  a  three  dimensional  search  for  q^,qg  and  <34 
was  performed  starting  with  qN'®  =  (1,0, 0,1,0). 


Table  2. 


N 

N 

<3l 

N 

<*3 

N 

<34 

4 

1.3876 

0.7441 

2.0007 

8 

1.4091 

0.9509 

2.0003 

16 

1.4124 

0.9774 

2.0003 

32 

1.4137 

0.9904 

2.0003 

True 

Values 

1.4140 

1.0 

2.0 
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